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Abstract 

Quantum Hall skyrmions are quantized solitons of a ferromagnetic 0(3) 
a— model. The reference, classical, solutions depend upon the interaction 
between the electrons and exhibit completely different asymptotic profiles for 
the physical Coulomb interaction than for the hard core interaction frequently 
used to generate variational wavefunctions. In this note we show, by means 
of numerical calculations on (large) finite size systems at v = 1, that this 
physically important difference, crucial for a sharp definition of their statistics, 
persists for the quantized skyrmions at v = 1. 
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Theoretical and experimental work over the last few years has adduced strong evidence 
that the physics of the quantum Hall (QH) system in the vicinity of (filling factor) v — 1 
is best understood in terms of a ferromagnetic ground state right at commensuration which 
supports charged excitations with non-trivial spin topology christened "skyrmions" 
The skyrmions arise as quantized solitons in a two dimensional 0(3) cr— model which de- 
scribes the low energy physics of the v = 1 state, and mutatis mutandis of other ferromagnetic 
QH states. Very general considerations of homotopy suggest the possibility of topologically 
stable and non-trivial classical solitons in any two dimensional 0(3) a— model 0] but the 
issue of their dynamical stability is more complicated. As was first noted by Belavin and 
Polyakov |J, the leading (gradient squared) term in their energy functional is scale invariant 
in two dimensions; hence the explicit solutions obtained by them exhibit an arbitrary scale 
parameter and algebraic asymptotic profiles. Consequently, dynamical stability, i.e. the 
existence of a well defined size for the skyrmions, depends crucially on the sub-leading terms 
which are system specific. 

In the theory of the u = 1 QH this issue can be resolved in two very different ways. 
For the physical, Coulomb, interaction, the effective a— model written in terms of the order 
parameter field n(r) is modified to [IJ 

= ^pi(n) • ftn - ^ s (Vn) 2 + ■ B - ^ / dV'|^ , (1) 

neglecting a Hopf term that does not affect the structure of the skyrmions. Here A is 
the vector potential of a unit magnetic monopole, p s the spin stiffness (p s = e 2 / (16^/2ne£) 
for v = 1), e the semiconductor dielectric constant, I the magnetic length and q(r) = 
n ■ (d x n x <9j,n)/47r is the topological density of the spin field. The Zeeman (third) and 
Coulomb (fourth) terms in £ e // break the spatial scale invariance of the gradient (second) 
term and their competition then sets the size and energy of the skyrmions which now depend 
on the dimensionless ratio g = (gpsB) / (e 2 / e£ 2 ) of the Zeeman energy to the Coulomb energy. 
On the other hand, one can study a "hard core" interaction, consisting of a delta function 
projected to the lowest Landau level, which is a standard device in QH theory to generate 
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variational wavefunctions |J. Remarkably, the effective a— model for this problem contains 
exactly the first three terms, i.e. it supports skyrmions of arbitrary size A at gusB = given 
by the Belavin/Polyakov solutions: 

n x = yfl - P(r) cos(^) , n y = ^1 - / 2 (r) sin(0) , n z = f(r) , (2) 

where f(r) = [(r/A) 2 — 4]/[(r/A) 2 + 4] @. For gusB ^ 0, the stable solution has vanishing 
size. Despite this last pathology of the hard core problem, it would appear that its solutions 
could again be used for a variational solution of the Coulomb problem, i.e. by optimizing A. 

It is a fact of some consequence, that this cannot be done. It turns out that the 
skyrmion profile in Eq (§) decays too slowly (n z — 1 ~ (A/r) 2 ) to yield a finite answer for 
the Zeeman energy, necessitating a calculation of its true asymptotic form in the Coulomb 
problem which can be shown to go as n z — 1 ~ e~ r ^, i.e. the Coulomb skyrmion has a 
much shorter tail than the hard core variant. Clearly, if this was not taken into account 
one would falsely conclude that Coulomb skyrmions are always microscopic. Further, it 
was shown by Yang and Sondhi |9j that the statistics of the skyrmions is ill-defined in 
hard core mult i- skyrmion states precisely on account of the long tails of the skyrmions but 
well defined for the shorter tailed Coulomb versions. Finally, the difference in asymptotic 
behavior was shown recently to lead to binding between singly charged skyrmions at small 



Zeeman energies, albeit at unphysically small values |10| . 

The above considerations were all in the framework of classical solutions of the a— model 
and the closely related Hartree-Fock approximation H. The corresponding quantum states 



require inclusion of fluctuations that quantize spin and angular momentum [fll|JT^1 and re- 
store the symmetries broken in the classical solutions, and one may wonder if the distinction 
between the tails of hard core and Coulomb skyrmions survives when this is done. In this 
note we report numerical calculations on the smallest skyrmions, i.e. those with one extra 
reversed spin for which the fluctuations are the most severe, and show that this important 
distinction is indeed preserved. 

Results: We study a QH system on a sphere || at one flux quantum less than the finite 



size version of v — 1, i.e. with flux quanta passing through the surface of the sphere 
and N — = 2K particles. This is the one quasielectron sector and its states are related 
to those in the one quasihole sector (N — + 2) by a particle-hole transformation. The 
spectrum in this sector has a set of states with equal values of the angular momentum L and 
spin S. The states with L = S = N/2 — l are the multiplet of the fully polarized quasielectron 
(it has just one reversed spin) and those with lower L, S are skyrmions of various sizes 
culminating in the infinite skyrmion (it has the size of the system) with L — S — 0. We 
have studied the smallest of the skyrmion states, i.e. with L = S = N/2 — 2 and maximum 
L z for simplicity for system sizes up to 300 electrons for both hard core and Coulomb 
interactions. 

For the hard core interaction, an exact wavefunction can be written down. In second 
quantized notation this takes the form, 

l*> = E{^=^ c ™tcli4- U + y^^ c ™t c Lu c kJh = i) , (3) 

where A = (2K + 1) £^= 2 1/n « (2K + l)(ln(2X) + 7 - 1), and 7 is the Euler's constant. 
From this we can easily compute the occupation of various orbitals by the down spins, 

n(K) = 1 - (1 - l/Nj/A 
n(K-l) = N < p(l-l/N <p )/A 

i (4) 
n(K -m) = (iV + l)/(mA) - (1 - l/N^/A m>2. (5) 

It is clear from these formulas that one of the down spins is localized at the north pole 
(top orbital). For N=300, in Fig 1 we plot the spatial density for the other down spin 
°"l( r ) — jpI2m=-K l^m(^) | 2n j.( m ) by excluding the m = K orbital. Here r = 2i?sin#/2 is 
the chord distance, R = \[Kl is the radius of the sphere. Evidently, cr^(r) is peaked near 
the nominal location of the quasielectron. In the same figure we plot the result for the exact 
wavefunction for Coulomb interactions which was obtained by numerical diagonalization. 
One can see a considerable difference in the degree of localization of the second reversed spin 
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between these two potentials. The difference in the tails, however, may not be immediately 
apparent from this plot. It becomes apparent when we plot the data on a double-logarithmic 
scale in Fig 2, where it is clearly seen that the down spins are strongly localized in the 
Coulomb case when compared with the hard core wavefunction; at the opposite pole there 
are seven orders of magnitude separating the two results(3.1 x 1CT 14 , vs. 8.5 x 1CT 6 )! 
Moreover, the hard core distribution follows the 1/r 2 decay expected from the sigma model 
analysis while the Coulomb result decays faster than any power consistent with the predicted 
exponential behavior. The asymptotic behavior on the sphere for the hard core potential 
to leading order in system size can easily shown to be ( j^ftp. ~ jc)i see the dashed 
line in Fig 1. Yet another probe of this difference is the large system behavior of the down 
spin occupation of a given orbital near the location of the quasielectron. The hard core 
wavefunction is easily seen to lead to a logarithmic vanishing of this quantity in the large 
system limit which implies that the quasielectron is not a well defined object at infinite 
system size: one of its reversed spins is spread over the entire system. By contrast the 
Coulomb quasielectron is well defined, examination of the variation of the occupation with 
system size displays clear evidence of finiteness in the large system limit (Fig 3). For 
comparison, we have produced a similar plot for the hard-core potential in Fig 4. Only 
the occupation of the top orbital shows scaling with system size consistent with its being 
localized. Others show noticeable deviations as expected for derealization of the second 
reversed spin. 

In summary, there is a crucial difference in the asymptotic behavior of hard core (scale 
invariant o— model) skyrmions and Coulombic skyrmions which is present also in the fully 
quantized and symmetric many-body states of the lowest Landau level problem. Hard core 
skyrmions have algebraic profiles and are ill-defined in the large system limit while Coulomb 
skyrmions have exponential profiles and have a definite size. This difference is crucial to 
the energetic stability of Coulomb skyrmions, their possessing a well defined statistics and 
to their binding at extremely small Zeeman energies. Two remarks are in order. First, the 
distinction we have drawn here between the skyrmions of the hard-core and the Coulomb 
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potentials trivially vanishes for infinite skyrmions as they have no tails. Second, any generic 
perturbation of the hard core interaction (one which produces a direct interaction in addition 
to the exchange) will break the scale invariance and give a well defined size to the skyrmions; 
the hard core limit, though important in that it enables an exact analysis, is non-generic in 
this regard. 
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FIGURES 

FIG. 1. The spatial distribution of the second reversed spin for the Coulomb and the hard core 
potential for a 300-electron size systems. The dashed line is the asymptotic behavior to leading 
order in system size given in text. The inset shows the tails of the distribution. 

FIG. 2. The log-log plot of Fig 1. The differences in the tail can be clearly seen. The density for 
the Coulomb tail is curved indicating it vanishes faster than any power. The straight line segment 
for hard core gives the 1/r 2 dependence. 

FIG. 3. The plot of occupation amplitudes for the top 4 orbitals for the Coulomb potential 
for up to N=300 plotted versus 1/N. All are essentially independent of system size and indicate 
localization of the second reversed spin. 

FIG. 4. Same as Fig 3 except for hard core potential plotted for comparison. Only the first 
reversed spin at the top orbital can be seen to be localized. The others show a noticeable decrease 
and ultimately vanish logarithmically with system size. 
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